// declare
void mfn(double A, double deltac, double q, double ftNL, char* name, float fNL);

double ratio3(double M, double A, double deltac, double q, double ftNL, double z);
double ratio4f(double M, double A, double deltac, double q, double ftNL, double z);
double ratio4(double M, double A, double deltac, double q, double ftNL, double z);
double ratio5(double M, double A, double deltac, double q, double ftNL, double z);
double ratio6(double M, double A, double deltac, double q, double ftNL, double z);
double ratio7(double M, double A, double deltac, double q, double ftNL, double z);

double nuc(double deltac, double M, double A, double q, double ftNL, double z);
double nucD(double deltac, double M, double A, double q, double ftNL, double z);

double M2R(double M, double A, double q, double ftNL);
double M3R(double M, double A, double q, double ftNL);
double M4R(double M, double A, double q, double ftNL);
double M5R(double M, double A, double q, double ftNL);
double M6R(double M, double A, double q, double ftNL);
double M7R(double M, double A, double q, double ftNL);

double M2RD(double M, double A, double q, double ftNL);
double M3RD(double M, double A, double q, double ftNL);
double M4RD(double M, double A, double q, double ftNL);
double d2R(double M, double A, double q, double ftNL);
double d3R(double M, double A, double q, double ftNL);
double d4R(double M, double A, double q, double ftNL);
double d2R1(double M, double A, double q, double ftNL);
double d2R2(double M, double A, double q, double ftNL);
double d3R1(double M, double A, double q, double ftNL);
double d3R2(double M, double A, double q, double ftNL);
double d4R1(double M, double A, double q, double ftNL);
double d4R2(double M, double A, double q, double ftNL);
double I2R1(double M);double I2R2(double M);double I3R1(double M);double I3R2(double M);
double I4R1(double M);double I4R2(double M);
double I2R1D(double M);double I2R2D(double M);
double I3R1D(double M);double I3R2D(double M);
double I4R1D(double M);double I4R2D(double M);
double He1(double x); double He2(double x);double He3(double x); 
double He4(double x);double He5(double x); double He6(double x);
double He7(double x);
double tgamma(double x);

double An(double n);
double Bn(double n);
/*
double hubble(double z);
double growth(double z);
double gr_int(double z);
*/
double const AC = 8.413e-10;
double const h = 0.7;
double const c = 299792.458;
double const Omm = 0.27;
